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Abstract 



In this paper we describe the asymptotic behavior, in the exponential time scale, 
of solutions to quasi-linear parabolic equations with a small parameter at the sec- 
ond order term and the long time behavior of corresponding diffusion processes. In 
particular, we discuss the exit problem and metastability for the processes corre- 
sponding to quasi-linear initial-boundary value problems. 
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1 Introduction 

ON ' 

Consider a dynamical system 

^ = W X% = xeR d , (1) 

together with its stochastic perturbations 

dX^ £ = b(Xt £ )dt + ea(X?' £ )dW u X^ £ = x G R d . (2) 

Here e > is a small parameter, W t is a Wiener process in M. d , and the coefficients a and b 
are assumed to be Lipschitz continuous. The diffusion matrix a(x) = (a,ij(x)) = cr(x)cr*(x) 
is assumed to be non-degenerate for all x. 

Let D be a bounded domain in M. d with infinitely smooth boundary dD. In this paper, 
with the exception of the last section, we assume that there is a point x G D such that 
for each x G D the trajectory of the dynamical system ([1]) starting at x is attracted to 
Xq. We assume that (b(x),n(x)) < for x G dD, where n(x) is the exterior normal to the 
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boundary of D. Let r e = min{t : X^' £ G dD} be the first time when X^' £ reaches the 
boundary of D. 

If e is small, then on any finite time interval the trajectories of the process X^ ,£ defined 
by (j2J) are close to the corresponding non-perturbed trajectory with probability close to 
one. Therefore, with high probability X^ ,£ enters a small neighborhood of the equilibrium 
point xq before leaving D. The process eventually exits D as a result of large deviations 
of Xt' £ from X% ([6], see also [8]). The large deviations are governed by the normalized 
action functional 

SoAv) = \ [ T Y, al3 (^t-HVt)m-b 3 (^))dt, T>0, peC([0,T],D), 
Jo i,j=i 

and the quasi-potential 

V(x ,x) = inf{<V(</>) = <P e C([0,T\,V),<p(0) = x ,<p(T) = x}, x G D. 

T,tp 

Here a u be the elements of the inverse matrix, that is a % i = (a -1 )^-, and So,t(¥>) = +00 
if if is not absolutely continuous. It is proved in [6] that e 2 In r £ converges in probability, 
as e I 0, to Vq = mm xe Q D V(x , x). Moreover, if the minimum V of V(x ,x) on dD is 
achieved at a unique point x* (which is true in the generic case), then X*f converges to 
x* in probability as £ J. 0. 

These statements imply various results for PDE's with a small parameter at the second 
order derivatives. In particular, consider the following initial-boundary value problem: 

i,j=l 1 1 

w e (0,x) = g(x), x G D, w £ (t,x) = g(x), t > 0, x G dD, (4) 

where g, for the sake of brevity, is assumed to be continuous on D. The case when 
u(t, x)\ xe QD = ip{x) with if> ^ g can be considered in a similar way. Assume that the 
minimum Vq of V(xq,x) on dD is achieved at a unique point x*. Let t : M + — > K. be a 
function such that t(e) x exp(A/e 2 ) as e [ with A > 0, that is hx{t{e)) ~ \/e 2 as s J, 0. 
Then 

\imw £ (t(e),x) = g{x ), if A < V , 

ej.0 

limw £ (t(e), x) = g(x*), if A > V , 

for x G D. 

Note that the solution to ([3])- (HI) can be expressed in terms of the transition semigroup 
associated with the family of processes Xf' £ , x G D. Namely, let Tfg{x) = Eg(X^ TS ), 
g G C(D). Then the function w e (t, x) = Tfg(x) is the solution to (jSI)-©- The semigroup 
T t e can be viewed as a small perturbation of the semigroup of shifts T t g(x) = g{Xf) 
associated with the dynamical system ([T]). 
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More general perturbations of T t may lead to nonlinear semigroups. Namely, consider 
the following problem: 

— = L £ u £ := e -Y] a l3 (x } u £ ) d ^ ^ + 6(s) ■ V x u £ (t,x) } x G D, t > 0, (5) 

at 2 ^ oxiOXi 

i,]=l ■> 

u £ (0,x) = g(x), x e D; u £ (t,x) = g(x), t>0, x G 3D. (6) 

When the coefficients are sufficiently smooth and the matrix a is positive-definite, the 
solution u £ exists and is unique in the appropriate function space (see Section [27TI below) . 
We can now define the semigroup Tl on C(D) via T £ g(x) = u £ (t,x), where u £ is the 
solution of ([S])-© with initial-boundary data g. 

For t > and x G D, we can define X*' x,e , s G [0, t], as the process which starts at x 
and solves 

dX*' x ' £ = b{Xl' x ' £ )ds + ea(Xl' x ' £ , u £ (t - s, X^ x ' £ ))dW s , s < t £ At, (7) 

r £ = min{s : X^' 6 G dD}, X'f' 6 = X% x ' £ , r e < s < t, 

where cr^, 1 < z,j 1 < d, are Lipschitz continuous and such that crcr* = a. The process 
X* ,a: ' £ will be called the nonlinear stochastic perturbation of ([T]). More precisely, Xl ,x ' £ 
corresponds to the nonlinear semigroup defined by As in the linear case, we have the 
following relation between u £ and the process Xl ,x,£ : 

u £ (t,x)=Eg(X t t ^). 

One of the important questions in the study of parabolic linear and nonlinear equations 
is the one concerning the behavior of solutions (or, in probabilistic terms, behavior of the 
corresponding diffusion process) as t — >• oo. In our case, when the small parameter e 2 is 
present in front of the second order term, the limit of u £ (t,x) as e — > 0, t — > oo, depends 
on the manner in which (e, t) approaches (0, oo). In the linear case this problem has been 
studied in [1] (see also [6], [8]). 

In Section [3] we study the asymptotic behavior of solutions to (jSJ)-® when e J. and 
t = t(e) x exp(A/e 2 ). As a first step, we shall introduce a family of linear problems 
which can be obtained from ([S])-® by replacing the second variable in the coefficients 
ciij in the right hand side of (|5J) by a constant c. The asymptotics of u £ can be then 
expressed in terms of the functions Vq(c) and g(x*(c)), where Vb(c) is the minimum of 
the quasi-potential of the linear problem and x*(c) is the point where this minimum is 
achieved. 

In Section H] we study the exit problem for the process Xl' x ' £ . We shall see that new 
effects appear when nonlinear stochastic perturbations are considered. In particular, even 
in the generic case, the distribution of the exit location X t '®£ need not be concentrated 
in one point. 

Some related problems concern the notion of metastability for nonlinear perturbations 
of dynamical systems with several equilibrium points. Let us consider the dynamical 
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system (pQ) in M. d and its perturbations As before, a and 6 are assumed to be Lipschitz 
continuous. Now we shall assume that the system has a finite number of asymptotically 
stable equilibrium points X\, Xk such that for almost every x G K d , with respect to the 
Lebesgue measure, the trajectory of (CQ) starting at x is attracted to one of the points 
x\,...Xk- We shall also assume that the vector field b satisfies (b(x),x) < A — B\x\ 2 
for some positive constants A and B. The case of more general asymptotically stable 
attractors (for instance, limit cycles) can be considered similarly, however for the sake of 
brevity we restrict ourselves to the case of equilibriums. 

The general theory of metastability was developed in [I] in the framework of large 
deviations (see also [6], [5], [8]). It was shown, in particular, that for a generic vector 
field h satisfying the assumptions above, for almost every x G M. d and A > 0, with 
probability which tends to one when e J, 0, the trajectory X^' E of (jSJ) spends most of 
the time in the time interval [0, exp(A/e 2 )] near a point x x G {x\, ...X}~}. This point is 
called the metastable state for the trajectory starting at x in the time scale exp(A/e 2 ). 
The metastable state can be determined by examining the values of the quasi-potential. 
Namely, let 

Vij = V{x u xj) = inf{So, T (^) = V e C([0, T], R d ), <p(Q) = x i: (p(T) = Xj }, 1 < i,j < m. 

T,ip 

These numbers determine a hierarchy of cycles along which the system switches from one 
metastable state to another with the growth of A ([!])• 

We can also study metastability for nonlinear perturbations of dynamical systems. It 
turns out that now the transition between the equilibrium points does not occur "imme- 
diately in the exponential time scale". This implies that now metastable states should 
be replaced by metastable distributions between the equilibriums. The description of 
metastable distributions is based on the study of the asymptotic behavior of solutions to 
flSJ)-© when e { and t = t(e) x exp(A/e 2 ). Note that metastable distributions also 
arise in pQ, [3], but for reasons which are different from what is discussed in this paper. 
Such a modification to the notion of metastability leads to a modified notion of stochastic 
resonance. 

We briefly address the problems of metastability in Section [51 where we also consider 
other generalizations and some examples. The issue of metastability in the case of an 
arbitrary number of equilibrium points and cycles will be addressed in a forthcoming 
paper. 

In this paper we considered nonlinear perturbations of a system with an asymptotically 
stable equilibrium. In this case, the exit from a domain containing this equilibrium occurs 
due to large deviations, and the exit time and exit distribution essentially depend on 
the perturbation. A related singular perturbation problem arises in the case when the 
equilbrium is stable but not asymptotically stable: for instance when the unperturbed 
system is Hamiltonian. Nonlinear stochastic perturbations in this case lead to a nonlinear 
version of the averaging principle. Say, in the case of one degree of freedom, the limiting 
slow motion is a diffusion process corresponding to a nonlinear operator on the graph 
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(compare with [6], Chapter 8) related to the Hamiltonian. We will consider these problems 
in one of the forthcoming papers. 

2 Preliminaries and Notations 

2.1 Quasi-Linear Equation 

Let D C R d be a bounded domain with infinitely smooth boundary dD. We shall say that 
/ : D — ► R belongs to C 2 (D) if / and all of its partial derivatives up to the second order are 
bounded and continuous in D. We shall say that a function / : (0, oo) xD->M belongs 
to C 1,2 ((0, oo) x D) if /, its partial derivative in t, and all of its partial derivatives up to 
the second order in x are bounded and continuous in (0, oo) x D. Note that a function 
/ G C 2 (D) can be extended to a continuous function on D and / G C 1,2 ((0, oo) x D) can 
be extended to a continuous function on [0, oo) x D. 

Let Oij = dji G C 2 (D x E), 1 < i,j < d, and b{ G C 2 (D), 1 < i < d. We also assume 
that there is a positive constant k such that k\£\ 2 < J2i j=i a ij( x i > x G D,u G R, 
£ G M d . Let g be an infinitely smooth function defined in a neighborhood of D. 

If e > and the coefficients a and b and the function g satisfy the assumptions 
listed above, then the equation (EJ)- (ED has a unique solution in the class of functions 
C 1 ' 2 ^, oo) x Dj_nC([0, oo) x D) (see Theorem 5, Chapter 6.2 of [7J). If g were to be only 
continuous on D, the existence and uniqueness of solutions to (jSJ)-© would hold in the 
class of functions which are locally C 1,2 -smooth inside (0, oo) x D and continuous up to 
the boundary. However, to simplify notations in later sections, we impose the smoothness 
condition on g. 

2.2 Action Functional 

Let a be a symmetric dx d matrix whose elements are bounded and Lipschitz contin- 
uous on R d and satisfy k\£\ 2 < J2i,j=i a ij( x )£,i£,j> x e £ e ^- d - Let a lj be the elements 
of the inverse matrix, that is a 1 ^ = (a -1 )^-, and a be a square matrix such that a = oo* . 
We choose o in such a way that are also bounded and Lipschitz continuous. 

Let Sq t be the normalized action functional for the family of processes X^' e satisfying 

dX*' e = b(X?> £ )dt + eo(X*> £ )dW t} 

where b is a bounded Lipschitz continuous vector field on R d . Thus 

SIM = \ f E oPfaM - UnWt - HVt))dt 

for absolutely continuous defined on [0, T], ip = x, and 5 , q T (<^) = oo if cp is not 
absolutely continuous or if <p ^ x (see [B]). Let V a (x,y) be the quasi-potential for the 
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family Xf' e in D, that is 

V a (x,y) = inf{^ T (^) : if G C([0,T],D),<p(Q) = x,<p(T) = y}, x,y ED. 

T,<p 

3 Asymptotics of the Solution 
3.1 Formulation of the Result 

Recall that (b(x),n(x)) < for x G dD, where n(x) is the exterior normal to the boundary 
of D. We shall assume that there is an equilibrium point Xq G D for the vector field b, 
and that all the trajectories of the dynamical system x(t) = b(x(t)) starting in D are 
attracted to xq. We also assume that there is r > such that (b(x), x — xo) < —c\x — xo\ 2 
for some positive constant c and all x in the r- neighborhood of xq. 

Let 5 > 0, D 5 = {x : x G D, dist(x,dD) > 5}, and u e be the solution of 
We shall be interested in the asymptotic behavior of u e (exp(X / e 2 ) , x), where A is fixed, 
x G D s , and e | 0. 

Let 

SW = mmg(x), g max = maxg(x), g 1 = min g(x), g 2 = maxg(x). 

x<=D x€D xedD x£dD 

Thus [5-1,5-2] C [5 min ,5max]- Let M : [5 min ,5 max ] -> M be defined by 

M(c) = min V a{ -' c) (x ,x), (8) 

xGdD 

where a(x, c) is extended to an arbitrary bounded Lipschitz continuous function satisfying 

^Kl 2 < Et=i Oii(x, c)Ui, Z e x G M d \ D. 

We next make some assumptions about the quasi-potential. It is not difficult to see 
that these assumptions are satisfied by a quasi-potential corresponding to generic a and b. 

We shall assume that for all but finitely many points c G [5 m in, 5max] the minimum 
in ([8]) is attained at a single point which will be denoted by x*(c). We assume that in the 
remaining points c 1 , c k the minimum is attained at two points of the boundary In this 
case the function x* : [g m i n , g ma , x ] — > dD is piece-wise continuous and has left and right 
limits at the points of discontinuity, as follows from the formula for the quasi-potential. 
Let x\{d) = lim cTcl x*(c) if d ^ g min and x* 2 (c l ) = lim cicl 2* (c) if d ^ 5 max , 1 < i < k. 
If c l = 5min, we define x\{c l ) as the point distinct from x^c 1 ) where the minimum of the 
quasi-potential is attained, and similarly we define x^c 1 ) if d = g max as the point distinct 
from x\(d) where the minimum of the quasi-potential is attained. 

We assume that x\(d) 7^ x^c 1 ), 1 < % < k (thus lim c j c i x*(c) 7^ lim^ x*(c) if d is an 
interior point of [g m m, 5max])- Define G\(d) = g(x*(c % )) and G 2 (d) = gi^x^d)). We can 
now define the piece-wise continuous function G : [5 m in,5max] — ► [51,52] via 

G{c)=g{x*(c)), ce[ grain ,g max ]\{c\...,c k }, G(c i ) = G 1 (d), 1 <i < k. 
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Let Co = g(xo) and define c\ as follows: 

If G(cq) > c , then c\ = inf{c : c > c , G(c) < c}. 

If G(co) < Co, then c\ = sup{c : c < Co, G(c) > c}. 
Note that c\ G [gi,^] since G(\g m i n , g max ]) Q [01,(72]- We shall require that the graph of 
G pass from the left of the diagonal to the right of the diagonal at C\. More precisely, we 
shall assume that if Ci > g mm , then for every 5 > there exists 5 G (0, S ] such that 

G(ci — 5)>cx — 5, 

and if c\ < g max , then for every <5 > there exists 5 G (0, <5o] such that 

G(ci + <5) < ci + 5. 

We also require that Co not coincide with any of the points of discontinuity c l for which 
Gi(c*) < < G 2 (c i ). 

Let A G (0, oo) and define function c(A) as follows: 

For < A < M(c ), let c(A) = c . 

For A > M(c ) and ci = c , let c(A) = c . 

For A > M(co) and C\ > Co, let c(A) = min{ci, min{c : c G [co, Ci], M(c) = A}}. 

For A > M(co) and Ci < Co, let c(A) = max{ci,max{c : c G [ci, Co], M(c) = A}}. 

Here we use the convention that the minimum of an empty set is +oo and the maximum 
of an empty set is — oo. (See Picture 1, where the thick line represents the graph of the 
function c(A) and A' is a point of discontinuity for the function c(A).) 

We also define A max = sup cg[co Ci] M(c) if c x > c and A max = sup ce[ci Co] M(c) if c x < c . 

Theorem 3.1. Let the above assumptions concerning the differential operator L e and the 
function G be satisfied. Suppose that the function c(A) is continuous at a point A G (0, oo). 
Then for every 5 > the following limit 

limM e (exp(A/e 2 ), x) = c(A) 
is uniform in x G D s , where u £ is the solution to |3P-([iJ). 

Remark. From Theorem 13. II and the definition of the function c(A) it follows that 

limM e (exp(A/£: 2 ), x) = C\ 

uniformly in x G D s if A > A max . Moreover, from the proof of Theorem 13.11 provided 
below it easily follows that the limit is uniform in (x, A) G D s x [A, oo) for each A > A max . 
Therefore, for each 5 > and A > A max there is Eo > such that 

\u e (t, x) — ci I < 5 

whenever e G (0,Eq), x G D 5 and t > exp(A/e: 2 ). 



7 




picture 1 

the thick line represents the graph of c(A) 



It is important to note that with the boundary values of g fixed, the limit C\ may 
still depend on the initial function through its value Cq at the equilibrium point. In the 
generic case, when the interval [g m in, <?max] can be represented as a finite union of intervals 
I\ U ... U J m , such that on the interior of each of the intervals the function G(c) — c is either 
strictly positive or strictly negative, the values of Cq belonging to the interior of the same 
interval will correspond to the same value of c\. 

The proof of Theorem I3.ll will use some properties of diffusion processes stated in the 
following section. Theorem I3.ll implies various results concerning the exit problem and 
metastability for the process X\' x,£ defined above. These questions will be considered in 
Sections @] and EJ 

3.2 Properties of the Diffusion Processes 

In this section we shall consider diffusion processes which are somewhat more general 
than those introduced in Section I2.2I Namely, we shall allow the diffusion matrix to be 
time dependent. The results stated in this section easily follow from the arguments of [5J, 
Chapter 4. 

Let a be a symmetric d x d matrix whose elements a^- are Lipschitz continuous on 
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R+ x M. d and satisfy 

d 

m 2 < ^^(i,x)e^<^iei 2 , (t,x)GR + xR d , eei d , (9) 

where A; and are positive constants. Let a be a square matrix such that a = a a* We 
choose a in such a way that are also bounded and Lipschitz continuous. 
Let Xt' e satisfy Xq' £ = x and 

dX?' £ = b{X^ £ )dt + ea(t, X?' £ )dW t , (10) 

where 6 is a bounded Lipschitz continuous vector field on IR d . Clearly, the law of this 
process depends on a only through a = aa*. 

For x G D, let t £ be the first time when the process reaches the boundary of D. Thus 
X*f is the location of the first exit of the process X*' 6 from the domain D. If a is close 
to a function which does not depend on time, then the asymptotics, as £ J. 0, of X*f and 
t £ can be described in terms of the quasi-potential. 

More precisely, let a be a bounded Lipschitz continuous matrix valued function on R d 
such that 

d 

m\ 2 < "ijWtej < ^iei 2 , x g R d , e g R d , 

where a = aa*. Let A be the set of points in dD at which mm xe QD V a (xo, x) is attained. 
This minimum will be denoted by v. 

Lemma 3.2. Suppose a is as above, and positive constants k and K are fixed. For every 
5 > there is positive x and a function p : M. + — > M. + with lim e |o p{s) = 0, such that for 
every a that is Lipschitz continuous, satisfies (TJJ) and 

SUp |cHy(t, x) — 0£ij(x)\ < x, (11) 

and every x G D s we have: 

(A) ' P(r e <exp((v + 5)/e 2 )) > l-p(e), 

(B) P(r £ > exp((v - 6)/e 2 )) > 1 - p(e), 

(C) P(dist(X r Y, A)<8)>1- p{e). 

The next lemma only requires the boundedness of the quadratic form a from above 
and below. 

Lemma 3.3. Suppose that positive constants k and K are fixed. There exists v o > such 
that for every < 5 < v o there is a function p : M + — ► R + with lim e |o p(z) — 0, such that 
for every a that is Lipschitz continuous and satisfies (TJJ) and every x G D 5 we have: 

(A) P(r £ >exp(V^ 2 ))>l-p(£), 

(B) P(|Xf ' £ -x \<5for all t G [exp^ 2 ), exp(^/5 2 )]) > 1 - p(e), 

(C) P(\X^ £ - x \ < 5 for allt e [0,exp(v /e 2 )}) >l-p(e). 
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An easy corollary of this lemma is that at an exponential time the process either can 
be found in a small neighborhood of Xq or has earlier crossed the boundary of the domain. 

Corollary 3.4. Suppose that positive constants k and K are fixed. For every 5 > there 
is a function p : M + — > M + with lim £ | p{t) — 0, such that for every a that is Lipschitz 
continuous and satisfies (TJ|) ; every x G D and t > exp(S/e 2 ) we have: 

P(|Xf' e -x \<5 or t £ < t) > 1 - p(e). 

Proof. Let 8\ > be sufficiently small so that there is a domain D with smooth boundary 
such that D Sl = D. If the process does not reach dD by the time t — exp(5/e 2 ), then we 
can apply Part (B) of Lemma 13731 to the domain D and the process starting at ^t-exp(s/e 2 )' 
and the result follows from the Markov property. □ 



3.3 Preliminary Lemmas 

The next step in the proof of Theorem 13.11 is to establish that u £ (exp(t / e 2 ) , x) is nearly 
constant on D 5 if t > is fixed and e is sufficiently small. This is accomplished in 
Lemma 13.71 below. 

Lemma 3.5. For every positive to and R there are positive C and Eq such that 

\u £ {t,x) -u £ (t,x )\ < Ce (12) 

whenever \x — Xq\ < Re, e < Eq and t > to- 
Proof Let v e (t,y) = u £ (t,x +sy), t G (0,oo), \y\ < 2R. Let Q = (0, oo) x B 2R {0) and 
Qo — (to, °°) x B R (0), where B r (0) is the ball of radius r centered at the origin. Then v £ 
satisfies the following partial differential equation: 

dv £ (t,y) lA. d 2 v £ (t,y) b(x + ey) ( . 

-^T~ = 2 ^ ^ y) ^dy~ + — ' VyV (t ' y) ' (t ' V) Q - 

Here Oy (t, y) = a-^XQ + ey, u £ {t, x + ey)) are uniformly bounded in e and satisfy k\^\ 2 < 
Y%j=ia*j(t,y)€i£j, t G (0, oo), y G B 2R (0), £ G R d . Moreover, sup y6B2fl(0) \b(x + ey)/e\ 
is bounded uniformly in e and |V y ay(t, y)\ can be estimated from above by a constant 
times 1 + | V y v £ {t, y) \ for (t, y) G Q, uniformly in e. Since the distance between Q and the 
boundary of Q is positive and \v £ \ is uniformly bounded in Q by max(|^i|, l^l); we can 
apply the a priori estimate (see Theorem 4, Chapter 5.2 of [7j or Theorem 6, Chapter 6.2 
of [7]) to bound sup( t) ^ e Q |V y f e | by a constant C independent of e. This immediately 
implies (112)1 . □ 
We shall need the following simple lemma about diffusion processes with the drift 
directed towards the origin. 
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Lemma 3.6. Let b be a C 2 smooth vector field on M. d such that (b(x),x) < —k\(x,x) for 
some positive k\ and all x G Mr. Let a(t,x) be a Lipschitz continuous function such that 
\(?ij(t, x)\ < k^, 1 < i, j < d, for all t > and x G M. d . Let Y x be the process starting at x 
that satisfies 

dY t x = b(Y x )dt + ea(t, Y t x )dW t . 

Then for every r,5 > there are R, sq and Eq, which depend on b and a only through k\, 
lt2, such that 

P(y s x \m s] efljte(O)) > 1-5 (13) 
holds for x G B r (0) and < e < e . 

Proof. Let h : R — > [0, 1] be a smooth even function with negative derivative on (1/2, 1), 
such that h(x) = 1 for < x < 1/2 and h(x) = for x > 1. Let 

f(t : x) = h(\x\ exp(— kit)/ Re). 

If R and s$ are sufficiently large and Sq is sufficiently small, then f(so\hie\,x) > 1 —5/2 
for x G B r (0), < £ < e . By the Ito formula, 

P(noim £ |e^ £ (o))>E/(o,F; o|ln£| ) = 

E/(s | \ne\,x) + E^ ° "'(£/ - ^)(s | Ine| - 8 , Y x )ds, 

where £ is the generator of the process Y t x . In order to estimate the integral in the right 
hand side, we note that 

(Cf-%)(t,x) > -C ma,x(h"(x))exp(-2k 1 t)/R 2 , 
at xe[o,i] 

where the constant C depends on k\ and k%. By taking R sufficiently large, we can bound 
the expectation of the integral from below by —5/2, thus proving (jTBl . □ 

Lemma 3.7. For every positive Xq and 5 there is positive Eo such that 

\u £ (exp(\/e 2 ),x) -u £ (exp(\/e 2 ),x )\ < 5 (14) 
whenever x G D s , e < eq and A > Ao- 

Proof. Choose r > small enough so that (b(x),x — Xq) < —c\x — xq\ 2 for some positive 
constant c and all x G B2 r {xo). First, let us prove a slightly weaker version of the lemma, 
namely that under the same assumptions (1141) holds for all x G B t (xq), e < Eq and A > Ao- 
Let X^' x ' e , s G [0, exp(A/e 2 )], be the process that starts at x and satisfies 

dx X,x,e = b( X ^ £ )ds+Ea(X*> x > £ ,u £ (exp(\/E 2 )-s,X^ x ' £ ))dW s , s < t £ Aexp(A/ £ 2 ), (15) 

t £ = mm{s : X*' x ' £ G dD}, X*' x ' £ = X^ x ' £ , r £ < s < exp(A/£ 2 ), 
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where e C 2 (D x 1), 1 <i,j < d, are Lipschitz continuous and such that oo* = a. (If 
the minimum in the definition of t £ is taken over an empty set, then it is considered to 
be equal to +00.) 

In order to avoid confusion, let us note that we switched the notation, using X^ ,x,£ 
with superscript A for the process on the interval [0, exp(A/£ 2 )] with exp(A/e 2 ) inside 
the coefficient in the right hand side of (|15p . The notation X l s ,x,£ with superscript t still 
applies to the process defined in ©. 

Given s > 0, we can define 

t[ = min(s |lne|,inf(s : \X^ X,£ - x \ = 2r)). 

From Lemma T3.3I applied to the domain .B 2 r(^o) and Lemma [3.61 it follows that for every 
5' > there are positive s , e and R such that 

P(X%" G B^xa)) > 1 - 6' (16) 

for all A > Ao, x G B r (xo) and e < sq. Using the Feynman-Kac formula, we can write 
M e (exp(A/e 2 ), x) as follows: 

u £ (exp(X/e 2 ),x) = Eu £ (exp(X/e 2 ) - t £ ,X^ x 

Therefore, 

\u £ (exp(\/e 2 ),x) -u £ (exp(\/e 2 ),x )\ <2 sup P(X^f' £ <£ B^xq)) sup \u £ (t,x)\ 

x£B r {x ) 1 t>0,xeD 

+ sup \u £ (exp(X/e 2 ) — So\\ne\,yi) — u £ (exp(X/e 2 ) — s \\ne\,y2)\- 

yi,y2&B Re {x ) 

The first term in the right hand side can be made smaller than 5/2 by ( |T6l) . while the 
second term can be made smaller than 8/2 by Lemma T3.5I Thus (fT4|) is true if we require 
that x G B t (xq). 

Next, let us prove the original statement. For v > 0, define 

r| = min(exp(f /e 2 ), r e ). 

From Lemma 13.31 it follows that there is < v < Ao such that 

1 (X^ £ G 5 r (* )) = 1, 

where the limit is uniform in A > Ao, x G D s . Using the Feynman-Kac formula, we can 
write M e (exp(A/e 2 ) , X^j clS follows: 

u £ (exp(X/e 2 ),x) = Eu £ (exp(X/e 2 ) - r|,X^f' e ). 

Therefore, 

\u £ (exp(X/e 2 ),x) - u £ (exp(X/s 2 ), x )\ < 2 sup P(X*f ' £ <£ B r (x )) sup \u £ (t,x)\ + 

xeD s 2 t>o,xeD 

sup |w e (exp(A/£ 2 ) - exp(f /e 2 ),yi) - u e (exp(A/e 2 ) - exp(v /e 2 ), y 2 )\. 
yi,y2€B r (x ) 

The first term in the right hand side tends to zero by ( JT71) . while the second term tends 
to zero since we proved that the lemma is true if we require that x G B r (xo). □ 
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3.4 Proof of the Theorem on the Asymptotics of the Solution 

In this section we prove Theorem I3.ll First, we examine the behavior of u £ for times 
which are small in the logarithmic scale. 

Lemma 3.8. There is a positive Vq such that for every < S < v there is Eo > such 
that 

\u £ (exp(X/e 2 ),x) - g(x )\ < 5 
whenever x G D s , < e < Eq and 5 < X < v . 

Proof. This lemma immediately follows from Lemma [3.31 □ 

The next three lemmas, central to the proof of Theorem 13.11 rule out certain types of 
behavior for the function u £ . 

Lemma 3.9. Suppose that ai < fi n < X n < a 2 for some constants ai,a 2 > 0, e n I as 
n — > oo, and 

u £n {exp(fi n /e 2 n ),x Q ) = /3i, u £n {exp(X n /e 2 n ),x ) = f3 2 
with (5\ 7^ /?2- Then there is 5 > such that 

exp(A„/4) - exp(fi n /e 2 n ) > ex^(5/e 2 n ) (18) 

for all large enough n. 

Proof. Consider the process X^ n,Xo,en given by fTTS]) with r £n being the first time when 
this process reaches the boundary of D. Define 

r 5 " = min(r £n ,exp(A n /4) - exp(/i n /e^)). 

Then 

u £ -(exp(X n /e 2 n ),x ) = Eu £ "(exp(X n /e 2 n ) - r £ ", X^ x °' £n ). 

The left hand side in this formula is equal to /?2- If ffTBT) does not hold, then by Part (C) of 
Lemma 13.31 and Lemma 13.71 the right hand side can be made arbitrarily close to (3\ along 
a subsequence. □ 

Lemma 3.10. Suppose that a± < fi n < X n < a 2 for some constants ai,a 2 > 0, e n j as 
n — > oo ; and 

u £n (exp(/j, n /e 2 n ),x ) = Pi, u £n (exp(X n /e 2 n ),x ) = (3 2 , (19) 

If Gram < A < Pi < ^max; then neither of the following is possible: 

(A) There is 5 > such that X n < M(/3 2 ) - 5, 

(B) There is 5 > such that G(c) < /3 2 - 5 for c e [/3 2 - 5, /3 2 + $] ■ 
If gmin < P2 < Pi < S'max; then neither of the following is possible: 

(A') There is 5 > such that X n < M(f3 2 ) - 5, 

(B) There is 5 > such that G(c) > (3 2 + 5 for c G \J3 2 - 5, fi 2 + 5] . 
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Proof. (A) Let us assume that A n < M((3 2 ) — 5. We can find A^ G [// n , A ra ] such that 
u £n (exp(\' n /e 2 n ),x ) = (3 2 and w(t,x ) < (3 2 for t G [exp(p n /e 2 n ),exp(X'Je 2 n )}. 

Let c%(x) = /3 2 ). Let r e be the first time when the process defined in ( TTOl) 

reaches the boundary of D. By part (B) of Lemma [3.2} we can choose x > such that 
whenever a is Lipschitz continuous and satisfies ([HD and (TTTT) . we have 

P(r £ " > exp(X'Je 2 n )) > P(r £ " > exp((M(/3 2 ) - S)/e 2 n )) > 1 - p(e n ), (20) 

for x G D s , where p does not depend on a and satisfies lim e j p(^) = 0. Choose > 
such that 

\aij(x,(3 2 ) - a,ij(x,P)\ < x 

whenever (3 G [/3 2 — 2(3' ', (3 2 + 2/3'], x G D 3 *. Choose a sequence p! n G [//„, A£J such that 
u £ -{p'Je 2 n ,x ) =p 2 -P' and w £ "(t,x ) G [/3 2 - /3 2 ] for t G [exp(/i;/4), exp(A;/4)]. By 
Lemma [3 .7\ we have 

\aij{x,(3 2 ) - aij(x,u £ "(t,x))\ < x (21) 

for x G -D^, t G [exp(/i^/e^), exp(A^/e^)], if e n is sufficiently small. Consider the process 
Xs n,X0,Sn given by (|15|) . with r £n now being the first time when this process reaches the 
boundary of D, and 

= min(r £ ", exp(X'Je 2 n ) - exp(p'Je 2 n )). (22) 

Then 

u^(exp(X'Je 2 n ),x ) = Eu £ "(exp(X>Je 2 n ) - f», x}k x °' £n ). (23) 

The left hand side in this formula is equal to /3 2 , while the right hand side can be made 
arbitrarily close to /3 2 — (3' by considering sufficiently small e n due to (l2"0l . Corollary 13.41 
(which applies due to Lemma [3. 9 p and Lemma [3.71 This leads to a contradiction. 

(B) Assume that G(c) < (3 2 — 5 for c G [/3 2 — 5, /3 2 + 5}. Let .4. be the set (consisting 
of either one or two points) where the minimum of V a ^''^ 2 \xo, x) is attained. From the 
definition of G it follows that g(A) C (— oo,/3 2 — 5). By part (C) of Lemma [3.2[ we can 
choose x > such that whenever a is Lipschitz continuous and satisfies © and (jTTj) . we 
have 

P(<?(X^' £ ) < & - 5/2) > 1 - p(e) (24) 

for x G -D 5 and all sufficiently small e. 

As in case (A), we can find /3' > and /i n < /i^ < A^ < A n such that 

u £n (exp(p'Je 2 n ),x ) = f3 2 - {3', u £ "(exp(X' n /e 2 n ), x ) = f3 2} 

and ([21]) holds for x G D H , t G [exp(/x^/e^), exp(A^/e^)], if e n is sufficiently small. We 
can again employ formula (1231 in which the left hand side is equal to f3 2 . The right hand 
side can be written as 

Eu £ ^exp(X'Je 2 n )-r £ ^xik X0 ' e n = 
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The first term in the right hand side here can be made arbitrarily close to P(r £n < 
r£n )(/?2 — (3') by Corollary 13.41 and Lemma [3.71 The second term on the right hand side 
can be estimated from above for large n by P(r £n = T £n )((3 2 — 5/4) due to (T24|) . This leads 
to a contradiction. 

The proof of (A') and (B') is completely similar to the proof of (A) and (B). □ 
Lemma 3.11. Suppose that A > 0. If c% > cq, then 

liminf(M £ (exp(A/e 2 ),x ) -c(A)) > 0. (25) 

ej.0 

If Ci < Cq, then 

limsup(M e (exp(A/e 2 ), x$) — c(A)) < 0. 

ej.0 

Proof. We shall only consider the first statement since the second one is completely sim- 
ilar. Note that 

liminf M £ (exp(A/e 2 ), x ) > c . (26) 

Indeed, otherwise by Lemma I3T81 there are (3 2 < Pi < c and sequences v < fi n < X n < A 
and e n I such that (fi~9j) holds. Note that the graph of G goes above the diagonal in 
a neighborhood of Cq, while (3\ and f3 2 can be taken arbitrarily close to Cq. Therefore, 
Part (B') of Lemma 13.101 leads to a contradiction. 

Thus, if ( 1251) does not hold, then is there are 5 > and a sequence e n { such that 

c - 5 < M £n (exp(A/e 2 ), x ) < c(A) - 5. 

We choose 5 sufficiently small so that the graph of G goes above the diagonal on the 
interval [cq — 25, c(A) — 5/2}. Take 5' > which will be specified later. 

For each c G [c — 5, c(A) — 5], by part (A) of Lemma l3~2l we can choose x(c) > such 
that whenever a is Lipschitz continuous and satisfies §§§ and (JTTJ) with aij(x) = aij(x,c) 
we have 

P(r £ < exp((M(c) + 5')/e 2 )) > 1 - p(e). (27) 
P(g(X:n > inf G(c) -5')>1- p(e) (28) 

c€[c-5' ,c+S'] 

For each c G [c — 5, c(A) — 5], find /(c) < 5' such that \ciij(x, c) — dij(x, c) \ < x(c) whenever 
c G[c- /(c), c + /(c)], x G D^ c \ 

Choose a finite subcovering of the interval [cq — 5, c(A) — 5} by the intervals (c m — 
l(c m )/2, c m +/(c m )/2) and take / = min m (/(c m )). Let c -5 = (3 < Pi < ••• < h = c(X)-5 
be such that A - A-i < V 10 > 1 < * < k - 

We claim that if i > 1, < A' < A and M £n (exp(A/e 2 i ), x ) G along a subse- 

quence, then M £n (exp(A'/e 2 ), x ) < for large enough n along the same subsequence. 
If this were not the case, then we would have 

M £ "(exp(A/e 2 ),x ) G and M £ "(exp(A7e 2 l ), x ) > A-i 
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along a further subsequence. Note that the function u £n (t,x ) must take values in the 
interval - V 10 > A + V 10 ] for * e [ ex P074)> ex P(V4)]> otherwise Part (B') of 

Lemma 13.101 leads to a contradiction. 

By the construction above and Lemma [3.71 there is m such that \c m — (3i\ < 5' and 

\aij(x,Cm) - aij(x,u £n (t,x))\ < x(c m ) 

for x G D^ Cm \ t G [exp(A'/e^), exp(A/e 2 )], if e n is sufficiently small. 

Consider the process X A ' x °' £n given by (1151) . with r £n being the first time when this 
process reaches the boundary of D, and 

r £ " = min(r £ ",exp(A/4) - exp(A;/4)). 

Then 

u £n (ex V (X/s 2 n ),x ) = Eu e "(exp(A/4) - r £ ",X^'°' £ "). 

The left hand side does not exceed If (exp(X/e 2 l ) — exp(A'/e^)) > exp((M(c m ) +5')/e 2 n ) 
(which is true if 5' is sufficiently small and n is sufficiently large), then from (1271) and ( 1281) 
it follows that the right hand side can be made larger than infc e [ Cm -<5',c TO +<5'] G(c) — 25'. This 
leads to a contradiction if 5' is small enough since G is a piece-wise continuous function 
which stays above the diagonal on [co — 25, c(A) — 8/2]. 

We have thus established that u £n (exp(A'/£: 2 ), 2o) < We can then extract a fur- 
ther subsequence such that u £n (exp(X / e^) , x ) belongs to one of the intervals [Pj-i, flj] 
with j < i. We can then take A" < A' and repeat the argument above to show that 
■u £n (exp(A"/e^), xq) < Pj-i- After at most k such steps, we obtain A < A such that 
u £n (exp(A/e 2 ), xq) < (3q along a subsequence, and A can be chosen to be arbitrarily close 
to A. This, however, is a contradiction with (1261) . □ 

Proof of Theorem \3.1\ By Lemma I3.7[ it is sufficient to prove that 

limw £ (exp(A/£ 2 ),x ) = c(A). (29) 

ej.0 

Case 1: < A < M(cq). Assume that 

limsupM £ (exp(A/e 2 ),xo) > cq. (30) 

ej.0 

Take Co < /3i < 0% < limsup £ i w e (exp(A/e 2 ), xq) such that M(/9 2 ) > A. By Lemma l3~8l 
there are sequences e n { and v o < \i n < A n < A such that (1T91) holds. Thus Part (A) of 
Lemma 13.101 leads to a contradiction with (1301) . The inequality 

liminf M £ (exp(A/e 2 ), xo) < Co 

eJ.0 

can be ruled out in the same way by referring to Part (A') of Lemma 13.101 
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Case 2: A > M(cq), c\ = cq. Assume that 

limsupw £ (exp(A/£: 2 ),Xo) > c . (31) 

ej.0 

Then, since G is piece-wise continuous and passes from the left of the diagonal to the 
right of the diagonal at ci, we can find 5 > and Pi, Pi such that 

c < Pi < (3 2 < limsupM £ (exp(A/e 2 ), x ) 

and G(c) < Pi — 5 for c E [Pi — 6, Pi + 5]. By Lemma 13. 8[ there are sequences e n [ 
and v < fi n < X n < X such that f|T9l) holds. Thus Part (B) of Lemma 13.101 leads to a 
contradiction with (l3T|) . The inequality 

liminf M £ (exp(A/e 2 ), xo) < Co 

can be ruled out in the same way by referring to Part (B') of Lemma 13.101 
Case 3: A > M(co), ci > cq. First assume that 

limsupM £ (exp(A/e 2 ), x ) > c(A). (32) 

We can repeat the arguments of Case 2 to show that ( 132]) implies that c(A) < c\. Then, 
since A is a point of continuity of c(A), we can find Pi, Pi such that 

c(A) < Pi < Pi < limsupti £ (exp(A/5 2 ), xo) 

and M(p 2 ) > A. By Lemma I3~8l there are sequences e n { and t>o < < A„ < A such 
that ( |T9|) holds. Thus Part (A) of Lemma [3.101 leads to a contradiction with ( 1321) . 
Finally, from Lemma 13.111 it follows that 

liminf M e (exp(A/e 2 ), x ) > c(A). 

Case 4: A > M(c ), Ci < c . This is completely similar to Case 3. □ 

Remark. If instead of the constant A in the argument of the function u £ in Theorem 13.11 
we have a positive function X(e) such that lim e | A(e) = A > 0, then 

limu e (exp(A(e)/e 2 ), x) = c(A). 
The proof of this statement requires only simple modifications to the proof of Theorem 13. 11 
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4 Exit From the Domain 



Let the differential operator U and the function G satisfy the assumptions of Section 13. 11 
Let x G D and A > 0. Recall that X^' x ' £ , s G [0, exp(A/e 2 )], is the process defined in ( fl5j) . 
with r e being the first time when this process reaches the boundary of D. We put t £ = oo 
on the event that the process does not reach the boundary by the time exp(A/e 2 ). Let 
T £ = min(r e , exp(A/e 2 )). Thus, if t £ < oo, then X^f' £ is the location where the process 
first exits the domain. Let p £ be the measure on D induced by X^f' e : 

p £ (A) = P(X^ £ G A), Ae B(D). (33) 

Let p £ be the restriction of p £ to dD: 

p £ (A) = P(X^ X ' £ G A), Ae 13 (dD). (34) 

Note that p £ is not a probability measure, since P(X*f' s G dD) < 1. In this section we 
shall examine the asymptotics of p £ and p £ when e J, 0. 

We shall distinguish several cases corresponding to different values of A. First consider 
the case when < A < M(c ). 

Lemma 4.1. If x E D and < A < M(c ), then lim e joP(^ e < exp(A/e 2 )) = 0. 
Proof. From Lemma 13.31 it follows that for 5 > there is < v o < A such that 

limP(X;°' x ' E G D for all < s < exp(v /e 2 )) = 1 (35) 

uniformly in x G D s . We claim that for each x > 0, 

\imu £ (t, x) = c uniformly in (t, x) G [exp(t> /e 2 ), exp(A/e 2 )] x D H . (36) 

Indeed, otherwise by Lemma [37fl we could find sequences e n j and A„ G (fo,A) such 
that either limsup^^ u £n (exp(X n /ef l ), xo) > cq or liminf^oo M £n (exp(A„/£: 2 ), xo) < Co. 
Suppose that the former is the case and that C\ = Cq (the argument in the cases when 
C\ > Cq and Ci < Cq is similar). Then, from the conditions imposed on the function G 
in Section |3~T1 it follows that there are Co < /3i < fa < hmsup^^ M e "(exp(A n /£ 2 ), x ) 
such that the graph of G goes below the diagonal in a neighborhood of the interval 
, /3 2 ] - Moreover, there are v < p' n < \' n < \ n such that u £n (exp(p' n / e^) , x ) = Pi and 
u £n (exp(X' n /ef l ), xo) = (32- This contradicts Part (B) of Lemma l3.10l thus establishing ( 1361) . 



From (1361) and Part (A) of Lemma I3T21 it follows that P(r £ < exp(A/e 2 ) — exp(t> /£: 2 )) 
tends to zero as e [ 0. From Corollary 13.41 it then follows that P(^exp(A/£ 2 )-cx P (?;o/£ 2 ) e 
tends to one as e | 0. By the Markov property of the process and due to fl35l) . this implies 
the statement of the lemma. □ 



Next, let us examine the case when A > A max . (Recall that A max = sup cg j C0 Cl ] M(c) if 
ci > c and A max = sup ce[ci co] M(c) if c x < c .) 
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Lemma 4.2. Suppose that x G D and X > A max . If the minimum of the quasi-potential 
mm x€ QD V a<y '' Cl \xo, x) is achieved at a single point x*{c\), then fx £ weakly converges to a 
probability measure \i concentrated at x*{c\). If the minimum is achieved at two points 
x*(c\) and x%{c\) and G\{c\) 7^ G^ci), then fi £ weakly converges to a probability measure 
fi concentrated at those two points. In this case fi(x*(ci))Gi(ci) + fi(x%(ci))G2(ci) = C\. 

Proof. Let A max < A' < A. Similarly to the proof of Lemma I4.1[ and using the fact that 
u e (exp(X/e 2 ), xq) converges to c\ for each A > A max , we can show that 

\imu £ (t,x) = ci uniformly in (t,x) G [exp(X' / e 2 ) , exp(A/e 2 )] x D*, (37) 

for each x > 0. Let A = {x*(ci)} if the minimum mm x€ g£) V a( -'' Cl \xo, x) is achieved at a 
single point, and A = {x\{c\) , x^Ci)} if the minimum is achieved at two points. Recall 
that t £ is the first time when the process X^ )X,£ reaches the boundary of D. From ( 1371) 
and Lemma 13.21 it follows that 

limP(r £ < exp(A/e 2 ) - exp(A7e 2 )) = 1 (38) 

and for every 5 > we have 

KmPCdistCX^^^^,^) < S) = 1. (39) 

This immediately implies the desired result for the case of a single minimum point. 

Let U 5 (y) C dD denote the S neighborhood of a point y on the boundary. In the case 
when the minimum is achieved at two points, we note that 

M £ (exp(A/e 2 ), x) = Ew e (exp(A/ e 2 ) — r e A (exp(A/ e 2 ) -exp(A'/e 2 )), X A f A ' ( £ cxp(A/£2) _ cxp{A , /£2)) ), 

where the left hand side tends to c 1; while the right hand side is equal to 

// £ (?7 5 (^( Cl ))) 5 (x*( Cl )) + if{U S {xl{ Cl )))g{x* 2 { Cl )) + a(s), (40) 

where \im e i a(e) = 0, as follows from (138!) and (!39l) . It also follows from (138]) and (|39l) 
that lim £l0 fi £ (dD \ {U s {x\{c x )) U U s (x*^)))) = 0, which, together with gOD, implies the 
desired result. □ 
Finally, we consider the case when c\ ^ cq and the function c(A) is continuous at a 
point A G (M(co), A max ). The cases c\ < Co and c\ > cq are completely similar to each 
other, so we shall only deal with the latter one. Let us introduce the needed notations. Fix 
Vq G (0, Co). By Theorem 13. 1[ lim e j w £ (exp(t> /£ 2 ), xq) = Co and lim e | w e (exp(A/e 2 ), xq) = 
c(A) > Co. For each c G [co, c(A)] we define 

A £ (c) = min(inf{A' > v : u £ {exp{X' /e 2 ),x ) > c},A). 



19 



Let a £ (c), c G [cq,c(A)], be the probability that the process X^' x ' £ reaches the boundary 
of D by the time exp(A/e 2 ) — exp(A £ (c)/e 2 ), that is 

a £ (c) = P(r £ < exp(A/£ 2 ) - exp(A £ (c)/£ 2 )). 

Since ot £ is left-continuous, it defines a measure v £ on B([cq,c(\)]) via z/ £ ([c, c(A)]) = 
a £ (c). It will be important to identify the limit of a 6 as e j 0. We define the function 
a:[c ,c(A)]^[0,l]by: 



a(c) = 1 — exp( 



C M dz 



z-G(z)' 



Since G is piece-wise continuous and its graph is above the diagonal in a neighborhood 
of [co, c(A)], the function a is a unique continuous function which satisfies the differential 
equation 

in the points of continuity of G and the terminal condition a(c(A)) = 0. Notice that 
a(c) G [0,1) for c G [co,c(A)]. The function a defines a measure v on i3([co, c(A)]) via 
u([c,c(\)]) = a(c). 

Lemma 4.3. 7/ci > c 0; c(A) is continuous at a point A G (M(c ), A max ) ; i/ien 

lima e (c) = a(c) 

for c G [c ,c(A)]. 

Proof. Assume first that the minimum of the quasi-potential is achieved at a unique point 
x*(c) for each c in a neighborhood of [co,c(A)], and therefore G is a continuous function 
there. Take 5,5' > which will be specified later and 5" > such that x, Xq G D s " . Let 
Po, Pi, Pk be such that Pq = Co, Pk = c(A) and < Pi — /3j_i < 5, 1 < i < k. 

Consider the processes Yp x ' E = X$ ^' x ' £ } 1 <i < k—1, and Y* ,x,£ = X^' x ' £ . Let r i,£ be 
the first time when the process Yg ,x ' £ reaches the boundary of D. Let B hX,£ , 1 < i < k — 1, 
be the event that r*' e < exp(A e (/3i)/e 2 ) — exp(A e (/5j_i)/e 2 ), and 5 fc,:c ' e be the event that 
r k > £ < exp(A/£ 2 ) - exp(A £ (/3 fc _ 1 )/e 2 ). 

Using Lemma 13.21 it is not difficult to show that for each 5' > there is 5 > such 
that 

limP^w n {dist(y r t/,x*(A)) > 5'}) = (42) 

uniformly in a; G -D 5 ". Since G is continuous, we can also make sure that 5 is small enough 
so that 

inn Pi /r "- n { //i v::;" ) - G(A)| > 5'}) = (43) 
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uniformly in x G D 5 " . We can write w e (exp(A/e 2 ), x) in two different ways 

u £ (exp(X/e 2 ),x) = 

E^(exp(A/, 2 ) — t £ A (exp(A/, 2 ) - e W (X £ (f3 t )/e 2 )), X^^^.^^) 

and 



u £ (exp(X/e ),x) = 

Eu £ (ex P (X/e 2 ) - r £ A (exp(A/, 2 ) - exp(A^(/3 J _ 1 )/e 2 )), ^(V^W.^)))- 
Upon subtracting the right hand sides of these two equalities, using the Markov property, 
Corollary 13 A\ Lemma [3.71 and (j43l . we obtain 

A(i-« £ (A)) = (a £ (A-i)-« £ (A))(G'(A) + ^i(^^)) + A-i(i-« £ (A-i)) + ^(^^), (44) 

where hi(i,e) < 5' and lim £ | /i 2 («, e) = 0. This implies the desired result once we recall 
that 5' and S can be taken arbitrarily small, since ( 144|) shows that a £ is a type of Euler's 
method approximation to the solution of ( )4Ti) . 

The condition of continuity of G can be easily removed once we recall that G may 
have at most finitely many points of discontinuity. □ 

Remark. Using similar arguments it is not difficult to show that 

limP(exp(A/e 2 ) - exp(A e (c )/e 2 ) < r £ < oo) = 0. 

Moreover, using (jHJ) and ( 1421) it is possible to show that in order to find the limit of p £ , 
one can take z/, which is the limit of v £ , and then take its push-forward by the function 
x* (since v is an absolutely continuous measure, it is not essential that x* may be unde- 
fined in a finite number of points). The push-forward of v will be denoted by p. Thus 
n(A) = u(c e [c , c(A)] : x*(c) G A), A G B(dD). 

Combining this with Lemmas 14.11 and 14.21 and Corollary 13.41 we can can formulate the 
following theorem. 

Theorem 4.4. Let p £ and \i £ be defined by [33^1 and p^ ), respectively. If x G D and 
< A < M(c ), then p £ — > 5 X0 , where 5 X0 is the probability measure concentrated at x 
and fi £ — > fx, where fi is the trivial measure, that is p{dD) = 0. 

If A > A max and the minimum of the quasi-potential mm x£ gjjV a( -' ' Ci \xq, x) is achieved 
at a single point x*(c\), then p £ and // weakly converge to a probability measure fi con- 
centrated at x*(ci). If the minimum is achieved at two points x\{c\) and x\(c\) and 
G\{c\) 7^ G2(ci) ; then p £ and p £ weakly converge to a probability measure p concentrated 
at those two points. Moreover, in this case p(xl(ci))Gi(ci) + /x(x2(ci))G 2 (ci) = c\. 

If Ci > Co and the function c(A) is continuous at a point A G (M(co), A max ), then take 
the measure v on B([cq, c(A)]) defined via z/([c, c(A)]) = a(c), where a is the solution of 
fi41\ ). The measures p £ weakly converge to the measure p which is the push-forward of 
v by the function x* . The measures p £ weakly converge to the measure c5 X0 + p, where 
c = l-p(dD). 
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Remark. This theorem still holds if instead of A in the definition of [i £ we have a positive 
function A(e) such that lim e j = A > 0. 

Corollary 4.5. If Ci > c and the function c(A) zs continuous at A G (M(c ), A max ), i/ien 
/or ei>en/ 5 > and x E D 5 we have 

limM e (exp(A/e 2 ),x) = / g(x* (c))dv{c) + g(x )(l - v([c , c(A)]), 

where v is the measure on £>([c , c(A)]) defined via z/([c, c(A)]) = a(c). 

Proof. The corollary immediately follows from Theorems 13.11 and 14.41 and the probabilistic 
representation of the solution to the initial-boundary value problem. □ 

5 Generalizations and Examples 

5.1 The Case of a Nonlinear First Order Term 

We could allow the coefficient at the first order term to depend on u e in ([5])-([6]): 

du £ (t, x) 



dt 



L £ u £ 



£ — Y] a i3 (x, u £ ) U y X) + (b(x, u £ ) + eh(x, u £ )) ■ V x u £ (t, x), x G D, t>0, (45) 

1,3=1 J 

u £ (0,x) = g(x), x G D, u £ (t,x) = g(x), t > 0, x G 3D. (46) 

All the assumptions made in Sections 12.11 and [3.11 remain in force, other than the following: 
instead of assuming that b is a vector valued function on D, we assume that b, b\ G 
C 2 (D x R), and there is a positive constant k' such that (b(x, u),n(x)) < —k' for x G dD, 
mGI, where n(x) is the exterior normal to the boundary of D. Moreover, we assume that 
for each u the vector field b(-,u) has a unique equilibrium point x which does not depend 
on u and that all the trajectories of the dynamical system x'{t) = b(x(t),u) starting in D 
are attracted to Xq. We now assume that there is a smooth function v defined on D, such 
that v(x ) = 0, v (x) > for x ^ x , and (b(x, u), Vv(x)) < —c\x — x \ 2 for some positive 
constant c, all u and all x. 

The definition of the function M(c) from Section [37X1 needs to be modified to allow for 
the dependence of the drift term on a parameter. Namely, now 

M(c) = min y a( '' c) ' b( ' c) (x ,x), 

xedD 

where 1/ a( ' c) ' 6( '' c) is the quasi-potential for the process whose generator is equal to 

— J2 aij (x, c) dx _ d ' + b(x, c) ■ V x u £ (t, x). 
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With this definition of M(c), Theorems 13.11 and 14.41 remain valid, and the proofs do not 
require serious modifications. 

5.2 Metastable Distributions in the Case of Two Equilibrium 
Points 

In this section we again consider the solutions w E to (E])-(l6l). Let all the assumptions about 
the domain D and the operator L e made in Sections 12.11 and [5.11 remain in force, except the 
following: instead of assuming the existence of a singe equilibrium point, we assume that 
there are two asymptotically stable equilibrium points Xi,x 2 E D such that for almost 
every x E D, with respect to the Lebesgue measure, the trajectory of (CQ) starting at x is 
attracted to either x\ or x 2 . Let D\ C D be the set of points in D which are attracted to 
x\ and D 2 C D the set of points attracted to x 2 . 

As before, we need to study the quasi-potential in order to determine the asymptotic 
behavior of u £ . While in the case of a single equilibrium, the function u £ was nearly 
constant in D 5 at times of order exp(A/e 2 ) (Lemma 13. T[) . now M £ (exp(A/e 2 ), x) will be 
close to M e (exp(A/e 2 ), xi) for x E D\ and close to M £ (exp(A/e 2 ), x 2 ) for x E D^. This 
explains why instead of freezing the second variable in the coefficients a%j in the right 
hand side of ([5]) , the way it was done is Section 12.11 now we put the variable equal to c\ 
in Di and c 2 in D 2 . More precisely, for Ci,c 2 E [^minj 5Wx]> let 

fcuc 2 (x) = c 1 xd 1 {x) + c 2 xd 2 (x) + (ci + c 2 )xd\(d 1 ud 2 )(x)/2, x E D, 

where xu is the indicator function of a set U C M. d . For a measurable positive-definite 
matrix- valued function a on D, we define 

V a (x,y) = M{3^ T {ip) : if E C([0,T],D),cp(0) = x,<p(T) = y}, x,y ED, 

T,tp 

where the normalized action functional S was defined in Section 12.21 Instead of function 
M used in Section [2~T| we now have functions M XljX2 , M X2>Xl , M Xu qd and M X2 ^r>. These 
are defined by 

M xl , X2 ( Cl ) = V a( -' f ^»(x h x 2 ), (47) 
M X2 , Xl (c 2 ) = V<>M-»(x 2 ,x 1 ), (48) 
M xljdD { Cl ,c 2 ) = min V<-' f ^^(x u x) t 

x£dD 

M X2 , dD (c 1} c 2 ) = min V^'^^ (x 2 , x). 

xedD 

It is not difficult to check that the right hand side of (1471) does not depend on c 2 and 
the right hand side of f|48|) does not depend on c\. For the process governed by equa- 
tion (jSJ), with dij(-, /ci,c 2 (')) instead of the nonlinear coefficients e ), the transition 
from x\ to a small neighborhood of x 2 occurs in time of order exp(M rijX2 (ci)/e 2 ) (pro- 
vided that the process does not exit the domain D earlier). Similarly, the transition from 
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x 2 to a neighborhood of x\ occurs in time of order exp(M X2tXl (c 2 )/e 2 ), while the transi- 
tion from X\ and x 2 to the boundary occurs in time of order exp(M E1) g£)(ci, c 2 )/e 2 ) and 
exp(M r2) a_ D (ci, c 2 )/e 2 ), respectively. 

In this section we would like to study the equation at a time scale which is sufficiently 
large for the process to make excursions between the neighborhoods of x\ and x 2 and 
back, yet not too large so that the process starting at x\ or x 2 does not exit the domain. 
Therefore we assume that 

max( sup M XuX2 (d), sup M X2:X1 (c 2 )) < M 9 , 

Cl G [flmiinSmax] C 2 £ [flmimSmax] 

where 

M 9 = min( inf M Ili9D (ci, c 2 ), inf M X2idD {c 1 ,c 2 )). 

Cl ,C2G[g m in,ffmax] Cl ,C 2 £ [SmimSmax] 

For example, if a and 6 are defined in the entire space M d , a is bounded and b satisfies 
(b(x),x) < A — B(x,x) for some positive constants A and B, then this condition will be 
satisfied for any domain D which contains a sufficiently large ball centered at the origin. 

Let ci = g(xi) and c 2 = g(x 2 ). Without loss of generality we may assume that c\ < c 2 . 
In order to formulate the results on the asymptotics of u £ , we need to introduce functions 
c 1 (X) and c 2 (A) which are similar to the function c(A) from in Section I3TT1 

Let < A < M 9 , and define c x (A) as follows: 

For < A < M xl)X2 {ci), let c x (A) = c x . 

For A > M XuX2 (ci), let c x (A) = min{c 2 ,min{c : c G [ci,c 2 ], M XuX2 (c) = A}}. 
Similarly, we define c 2 (A) as follows: 

For < A < M X2tXl (c 2 ), let c 2 (A) = c 2 . 

For A > M r2j:Cl (c2), let c 2 (A) = max{ci,max{c : c G [ci, c 2 ], M X2iXl (c) = A}}. 

Let A* = inf{A : c x (A) > c 2 (A)}. Assume that at least one of the functions c 1 and c 2 
is continuous at A*. Let c* = c l (X*) if c 1 is continuous at A* and c* = c 2 (A*) otherwise. 
Let c x (A) = min(c 1 (A), c*) and c 2 (A) = max(c 2 (A), c*). 

We can now formulate the following analog of Theorem 13.11 

Theorem 5.1. Let the above assumptions be satisfied. Suppose that the function c l (X) is 
continuous at a point A G (0, M 9 ). Then for every 5 > the following limit 

limu £ (exp(X/e 2 ),x) = c 1 (A) 

is uniform inx G Df . Suppose that the function c 2 (A) is continuous at a point A G (0, M 9 ) . 
Then for every 5 > the following limit 

limM £ (exp(A/e 2 ),x) = c 2 (A) 

is uniform in x G D%. 
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Remark. If A > A*, then c 1 (X) = c 2 (A) = c*. It is not difficult to see that the limit 

liniM e (exp(A/e 2 ), x) = c* 

e|0 

is uniform in (x, A) G D s x [A, oo) for each A > A*. Therefore, for each 5 > and A > A* 
there is e > such that 

\u £ (t,x) -c*\<5 

whenever e G (0, £o), x G D 5 and t > exp(A/£ 2 ). 

Recall that X^ £ , s G [0,exp(A/e 2 )], is the process defined in (fl~5l) . t £ is the first 
time when this process reaches the boundary of D and t £ = min(r £ , exp(A/£ 2 )). Since we 
assume that A < M 9 , the probability that r £ < exp(A/e 2 ) now tends to zero as £ J, 0. The 
distribution of the random variable X^f ,£ will be concentrated near the points X\ and x 2 . 

Theorem 5.2. Suppose that c\ ^ c 2 . If the function c^A) is continuous at a point 
A G (0, M 9 ) and x G D\, then the distribution of the random variable X^f' e converges to 
the measure /i^ = a\8 Xl + a 2 5 X2 , where the coefficients a\ and a 2 can be found from the 
equations c x (A) = a\C\ + a 2 c 2 , Oi + a 2 = 1. 

// the function c 2 (A) is continuous at a point A G (0, M 9 ) and x G D 2 , then the 
distribution of the random variable X r f' £ converges to the measure fi 2 = a\5 xi + a 2 8 X2 , 
where the coefficients a\ and a 2 can be found from the equations c 2 (A) = a\C\ + a 2 c 2 , 
at + a 2 = 1. 

If A G (\*,M 9 ) and x G D, then the distribution of the random variable X^f' £ con- 
verges to the measure //* = ai5 Xl + a 2 5 X2 , where the coefficients a\ and a 2 can be found 
from the equations c* = a\C\ + a 2 c 2 , a% + a 2 = 1. 

The proofs of Theorems 15.11 and 15.21 rely on the same techniques as those used in the 
proofs of Theorems 13.11 and I4.4[ and therefore will not be presented here. If a is bounded 
and b satisfies (b(x),x) < A — B(x,x) for some positive constants A and B, then similar 
results can be formulated for the Cauchy problem in M. d and the corresponding nonlinear 
perturbations of the dynamical system in ~R d . In this case we do not need the condition 
A < M 9 , but can instead consider all A G (0, oo). 

5.3 Examples 

In this section we give two examples when we can easily calculate the function M(c) 
defined in Section I3TT1 

In the first example, we assume that the domain D is one dimensional: D = (A, B). We 
assume that a(x, u) G C 2 ([A, B] x E) is positive, b(x,u) G C 2 ([A,B] x R), db(x,u)/dx < 
k < and there is a point xo G (A, B) such that b(xa,u) = for all u. In this case the 
operator 

L £ u = — a(x, u)u" + b(x, u)u' 
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satisfies the assumptions of Section 15.11 The quasi-potential, which will now be denoted 
by V c , is given by 

V°(xo,x) = -2 f X ^\dy, 
Jx a{y,c) 

as is easily seen from the definition of the action functional (see Section [272T) . Therefore, 
M(c) = min(^o, A), V c (x , B)) = min(-2 t ^\dy, -2 f ^\dy). 

The function G(c) may take at most two values: g(A) and g{B). In particular, the value 
g(A) is taken on the set {c : V c (xq, A) < V c (xq, B)}, while the value g{B) is taken on the 
set {c : V c (x ,A) > V c (x ,B)}. 

Picture 2 shows an example of the graphs of functions M(c) and G(c) in the case when 
g(A) = ini X £[A,B] g(x) < sup^g^ ^j g(x) = g(B). From Theorem 14.41 and the discussion in 

Section 15.11 it follows that for A > A max , the distribution of the random variable X^f' e 
converges to the probability measure /i concentrated at the end points of the segment. 
This measure can be found from the relations 

(i(A)g(A) + (i(B)g(B) = c h n{A) + fi(B) = 1. 




26 



In the second example we assume that D C M. d contains the origin xq = 0. Let the 
operator L £ be as follows 

i£ « = y£ Gy ' (M) Slr + {A{u)x) ' Vm ' 

i,j=l 1 ] 

where a is a positive-definite matrix which depends smoothly on u, and A is a matrix with 
negative eigenvalues which depends smoothly on u. As has been demonstrated in [2], the 
quasi-potential, which we shall denote by V c , is given by the quadratic form 

V c (x ,x) = )-{B~ l (c)x,x), 
where the matrix B is given by 

P CO 

B(c)= I exp(A(c)t)a(c)exp(A*(c)t)dt. 
Jo 

Therefore, M(c) = mm x£ Q£,(B~ 1 (c)x,x)/2. 
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